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Abstract 

One-loop radiative corrections to the leptonic tensor in high energy bremsstrahlung 
on heavy nuclei are calculated. Virtual and real photon radiation is taken into account. 
Double bremsstrahlung is simulated by means of Monte Carlo. Numerical results are pre- 
sented for the case of muon bremsstrahlung in conditions of the COMPASS experiment 
at CERN. 

1 Introduction 

Charged lepton Bremsstrahlung on nuclei has been studied both theoretically and experimen- 
tally for many years (see e.g. textbooks [U12] and references therein). This process contributes 
to energy losses of a lepton propagation through matter, which is relevant for many applica- 
tions. 

The similar process of high energy pion bremsstrahlung is used for extraction of pion 
polarizability [Hill]. In the modern COMPASS experiment j5j [6] the muon bremsstrahlung is 
used as a reference cross section and for estimates of systematic uncertainties. For this reason 
differential distribution of the muon bremsstrahlung should be predicted with high accuracy. 
That requires to take into account several effects beyond the tree-level Born approximation. 
So one needs to consider multiple photon exchange with the nucleus (important for large Z 
values), electromagnetic nuclear elastic and inelastic form factors, screening of the nucleus 
by the electrons surrounding it, and inelastic interactions of the projectile particle with the 
atomic electrons (see Ref. [8] and references therein). Besides those, we have to take into 
account also the vacuum polarization in the exchanged photon and at least one-loop radiative 
corrections to the lepton tensor. In this paper a new calculation of the latter is presented. 

Ref. [H] gives a comprehensive report on the calculation of one-loop corrections to virtual 
Compton scattering (ep — > epj). The lepton bremsstrahlung on a heavy nucleus we met here 
is a specific case of the general problem. Results of Ref. [9] for the Compton tensor also 
can't be directly applied to the problem under consideration, since the actual kinematical 
conditions (see Eq.Q below) deserve a special treatment. The latter includes keeping an 
exact dependence on the lepton mass and providing numerical stability of the corresponding 
computer code. 
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Since for this kinematical region we have to keep the exact dependence on the lepton mass, 
analytic formulae for the corrections become lengthy and cumbersome. Moreover to provide 
a possibility to impose various experimental cuts, we perform integration over the final state 
phase space numerically. 

The paper is organized as follows. In the next section we give the notation and the explicit 
expression for the Born cross section. Sect. [3] presents the calculation of various radiative 
correction contributions: the one due to a single virtual loop, the one due to additional soft 
photon emission, and the one due to double bremsstrahlung. Numerical results and conclusions 
are given in the last section. 



2 Preliminaries 

At the Born level we can represent the differential spectrum of the hard photon produced in 
the process 

Z(pO + A(P) -> /(pj,) + 1 (k) + A(P') (1) 

in the form [1, ,2] 
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where ui = k° is the emitted photon energy; ip is the azimuthal angle of the scattered lepton; 
Z is the nucleus charge; m is the lepton mass; £1(2) and Pi(2) are the energies and 4-momenta 
of the projectile (scattered) leptons, 

ci, 2 = cos(fcp 12 ), si )2 = sin(fcp 12 ), xi,2 = kpi, 2 , Q 2 = -(pi - P2 - k) 2 . (3) 

Here and in what follows, it is assumed that the lepton mass is small compared with the atom 
mass, while the energies are large: 

m < M A , Ex > m, E 2 > m, > m. (4) 

Let us rewrite the Born cross section (T2J) via a set of form factors: 



dcu 27T 



J dci y dc 2 y d ^ f^i 05 ^* ~ ^u(Qs,t s ,u s )Ei 



-1 -1 

- ^22 ( ( ls,ts,u s )E 2 + J r [ 2 \q s ,t s ,u s )E 2 E 1 + J r !£\q s ,t s ,u s )E 1 E 2 ^j, (5) 

<?s = 2x2 - 2xi + Q 2 , t, = 2xi - m 2 , m s = -2x 2 - m 2 , 

where e is the electron charge. The notation in the above expression is adjusted to the one used 
in the SANC [11] system, where the relevant expressions can be found as for the Born-level 
form factors as well as for the ones in the one-loop approximation. 
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Studying the differential distribution in the scattering angles of the Born cross section, one 
can see that it is peaked in the kinematical domain, where 



• — - _ m 

«Pl,2 ~ PlP2 ~ TT- ( 6 ) 

•C'l 

For the case of high energy muon scattering {E\ ~ 100 GeV) being under consideration now, 
the angles become small. Moreover, one should be careful with the dependence on the lepton 
mass, since m 2 ~ xi,2 i n this domain. On the other hand, we can safely drop some terms, 
proportional to the small ratio m 2 /Ef. As concerning the momentum transferred, contrary to 
the case of the Rutherford scattering, it can't go down below the kinematical threshold value 



,2 



1 ~ Tfl LO 

Qmin = Y Q m in = 2E1E2 ' 

In the ultra-relativistic approximation applicable in our case, after an integration over the 
whole phase space, one gets the Born-level photon spectrum in the simple form 

E\ + E 2 — — E1E2 J , to = E\ — Ei- (8) 



du mfaEfy 1 z 3 



3 One-Loop Corrections 

We subdivide the contributions of the one-loop QED corrections into three parts: 1) the one 
due to a single virtual loop; 2) the one due to soft real photon emission; 3) and the one due 
to additional hard photon emission (double bremsstrahlung) . 



3.1 Virtual Loop Contribution 

Representatives of the Feynman diagrams corresponding to the first type of corrections are 
shown in Fig. 13.11 This contribution was computed with help of the automatized computer 
system SANC [11] . The system provided the set of form factors calculated keeping the exact 
dependence on the lepton mass. The form factors are expressed via a number of one-loop 
master integrals (Passarino-Veltman functions), which are called from a SANC library. The 
infrared divergence in the relevant integrals is regularized by a fictitious photon mass A. So, 
the virtual loop contribution takes the form 

2 3 1 1 2n 

-1 -1 

^(q s , t s , u s )El + F${q„ t s , u s )E 2 E 1 + T$ \q a , t s , u s )E l E 2 ^j . (9) 

Here we adopt the SANC notation for the arguments and the normalization of the form factors. 
More details about the evaluation of form factors for such processes within SANC can be found 
in Ref. pj]. 
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Figure 1: Representatives of Feynman amplitudes with single virtual loops. 



3.2 Soft Photon Contribution 



Using the phase space slicing method we define the soft photon contribution as the one of 
the process with emission of an additional photon with energy below a certain value Co, which 
is small compared with the beam energy In our calculations we perform the spicing in the 
laboratory reference frame where the nucleus is at rest. Using the standard techniques of soft 
photon emission calculations we get the corresponding correction, which is factorized before 
the Born cross section: 
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(10) 



The infrared divergence of the soft photon contribution is regularized by means of a fictitious 
photon mass A, the same as in the virtual loop contribution. One of the internal cross checks 
of the calculation is the cancellation of the dependence on this auxiliary parameter in the the 
sum of two contributions. 



3.3 Double Bremsstrahlung Contribution 

Here we start with the completely differential expression for the matrix element squared. 
Some of the Feynman diagrams for this process are shown in Fig. 13.31 The two photons are 
treated in a symmetric way. In particular, the condition ui i2 > Co is applied for both the 
photons. The identity factor 1/2! is taken into account. Cancellation of the dependence on 
the parameter Co is checked numerically in the sum of the soft and hard contributions. The 
contributions of double real photon emission is computed by means a Monte Carlo integrator 
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based on the VEGAS algorithm [13] . Seven- fold integration over the whole final state phase 
space (including integration of the detected photon energy) is performed. The distribution in 
the detected photon energy is extracted in course of the integration using weights provided by 
VEGAS for each thrown kinematical point. 




Figure 2: Representatives of Feynman amplitudes for double bremsstrahlung. 



4 Numerical Results and Conclusions 

Summing up the considered above contributions we get the 1-loop corrected cross section in 
the form 

da Corr ^ d(J Born ^ d(J Virt ^ ^Soft ^ da Hard 

duj du dui du duj 

In Table H] there are numerical results for the specific contributions obtained for the following 
set of conditions: 



E 1 = 190 GeV, Z = 82, Q max = 0.0075 GeV 2 , 

u) = 0.001 GeV, M max = 3.75 -m^, P± n = 0.045 GeV, 

mi = m M = 0.10566 GeV, (12) 



where M max is the maximal allowed invariant mass of the muon plus hard photon final state 
system; P^ in is the minimal allowed transverse momentum of the outgoing muon; Z is the Pb 
nucleus charge. For the sake of simplicity, while computing the numbers for the Table we put 
a simple cut on the second hard photon energy: U2 < 0J\. By subscripts 1 and 2 we denote 
the results obtained with Co = 10~ 3 and uj = 10~ 4 , respectively. The relative corrections 5\ t i 
are computed as 

Si 2 = : n ' ■ 100%. (13) 

For a realistic simulation of spacial resolution and cluster energy threshold of the COM- 
PASS calorimeter in addition to the conditions (j!2p . we apply the following treatment of events 
with two hard photons: 

1) max(a>i, uj 2 ) > i.e. at least one of the photons should have an energy exceeding the 
threshold; 

2) if both the photons have energies above the threshold and the angle between their 
momenta is more than # 77 , the event is dropped; 
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u/E, 


Born 


Virtual 


Softx 


Hardi 


*L, % 


Soft 2 


Hard2 


s 2 , % 


0.3 


15677(1) 


76.8(4) 


- 260.1(1) 


226.9(3) 


+0.28 


-307.0(1) 


273.7(3) 


+0.28 


0.5 


10836(1) 


77.9(2) 


- 319.0(1) 


280.0(3) 


+0.36 


-377.4(1) 


338.1(3) 


+0.36 


0.7 


7337.7(1) 


76.9(2) 


- 363.3(1) 


297.1(2) 


+0.15 


-430.9(1) 


364.8(2) 


+0.15 


0.9 


1267.4(1) 


20.5(1) 


- 111.1(2) 


65.9(1) 


-1.95 


-132.4(2) 


87.2(1) 


-1.95 



Table 1: Contributions to the corrected differential cross section in pbarn/GeV versus the 
photon energy fraction. 



3) if the angle between their momenta is less than # 77 , the reconstructed photon energy is 
taken as the sum of the two: u = Ui + u 2 ; 

4) if one of the photon energies is below the threshold, the reconstructed photon energy is 
taken as the sum of the two: uj = uj\ + uj 2 . 

The parameter values correspond to one of data analysis procedures used by the COMPASS 
experiment, 

c^th = 7GeV, 6+y = 3 mrad. (14) 
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Figure 3: Relative contribution of one-loop corrections for realistic set-up vs. the photon 
energy fraction. 

For the realistic set-up, the size of the resulting correction is found to be below the one 
percent level. That is due to the fact that the correction is proportional to a/ (2ir), and in our 
case there is no any enhancement factors. In particular, even so that the beam energy is so 
large compared with the muon mass, the contributions of the order O [aln(Ef / mt)) cancel 
out in the sum of different contributions due to destructive interference of the initial and final 
state radiation. As can be seen from the Table H] at the end of the spectrum (uj — > Ex), where 
the phase space of additional hard photon emission is vanishing, we have a negative peak of 
the resulting radiation correction, which behaves there like aln((E% — u)/E\). But this peak 
is effectively washed out from the end of the spectrum in Fig. HI because of the the additional 
conditions ffT4"j) on event selections. 
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An analogous study was performed for the case of pion bremsstrahlung in Ref . [H] , where 
a similar behavior and magnitude of the one-loop corrections have been obtained within the 
scalar QED. 
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